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ABSTRACT 

Atmospheric circulation on tidally-locked exoplanets is driven by the absorp- 
tion and reradiation of heat from the host star. They are natural heat engines, 
converting heat into mechanical energy. A steady state is possible only if there 
is a mechanism to dissipate mechanical energy, or if the redistribution of heat 
is so effective that the Carnot efficiency is driven to zero. Simulations based 
on primitive, equivalent-barotropic, or shallow-water equations without explicit 
provision for dissipation of kinetic energy and for recovery of that energy as heat, 
violate energy conservation. More seriously perhaps, neglect of physical sources 
of drag may overestimate wind speeds and rates of advection of heat from the 
day to the night side. 

Subject headings: hydrodynamics — waves — binaries: close — planetary systems — 
stars: oscillations, rotation 



Introduction 



Showman et al.l (120071 ) have called attention to the divergence in results among at- 
tempts to simulate atmospheric ciculation on strongly irradiated jovian exoplanets ("hot 
Jupiters"). Simulations that allow explicitly for absorption and rejection of heat by the atmo- 
sphere, whether by simple Newtonian-cooling schemes or by radiative transfer, generall y find 



prograde equatorial wind speeds 1-5 km s ^ (IShowman fc Guillot 



Cooper fc Showmanll2005l : iLangton fc Laughlinll2008l : IShowman et al. 



2002 



Burkert et al. 



2008a 



20051: 



Dobbs-Dixon fc Lin 



20081 ). For comparison, the sound speed of nondegenerate molecular hydrogen is approxi- 
mately 2.4(r/10^K)i/2 kms'S so many of these winds are supersonic. Simulations that rep- 
resent uneven heating indirectly, via prescribed lateral vari ations in the thickness or depth of 
the c i rculating layer, ty pically have peak winds < 1 km s~^ (ICho et al.ll2003l : iLangton fc Laughlin 
20071 : ICho et al.l 120081). and they some times predict th a t the equatorial flow is retrograde. 
In the simulations of ICho et al.l (120031 ) and ICho et al. J2008h . the root-mean-square wind 
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speed is put in by hand in the initial conditions. Unhke wind speeds, hght curves of 
these planets have been directly measured, or at l east constrained, by infrared observations 
( iHarrington et al.ll2006l . 120071 : ICowan et al.l 120071 : iKnutson et al.l 120071 ). The models gener- 
ally do not predict these light curves very well: in particular, they sometimes produce a 
60-90° prograde shift of the photospheric temperature extrema with respect to the substellar 
and antistellar points (ICooper fc Showmanll2005l : iLangton fc Laughlinll2007l ), and they some 



times predict shifts of the wrong sign, as for the cold spot on HD 187933b (IShowman et al. 
2008ah . 



There may be many reasons for the disagreements among the models an d between the 



models and th e observations. For e xampl e, with the exception of those by iBurkert et al. 



(120051 ) and by iDobbs-Dixon fc LinI (120081 ) . the calculations use approximations based in 
part on the presumption of subsonic flow, even though the results sometimes violate that 
presumption. The purpose of the present paper, however, is to emphasize one aspect of 
the physics that has not been adequately respected by any of the simulations: namely, that 
atmospheric circulation is a natural heat engine, and that dissipation of mechanical energy 
must therefore be explicitly addressed. While it may not be easy to decide how dissipation 
occurs on real planets, it is possible that if a common model for dissipation were adopted, 
then the calculations might agree better on gross features of the flow such as peak wind 
speeds and light curves. Without such a common model, convergence to a common result 
seems unlikely. 



Production of mechanical work 



Three of the most recent models of radiatively driven circulation on exoplanets have been 
based on d irect integration of the so-called "primitive equatio r is" adapted from terrest rial me- 



teorology (IShowman et al.l l2008al : iMenou fc Rauscher 



20081 : IShowman et al. 2008b). Com- 



pared t o previous efforts to model the dynamics of atmospheric circulation, IShowman et al. 
(l2008bl ) in particular have significantly improved the model for coupling the flow to the ra- 
diation field. Our discussion does not do justice to those improvements, because our concern 
is with the dynamical part of the model, which has much in common with other efforts, 
including some of those that employ height-averaged versions of the primitive equations: 
more detail is given in §3 



The primitive equations are 



Dt 



-Vh$ - fkx Vh, 
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dp 



1 

P 
duj 
dp 



0, 



(2) 
(3) 
(4) 



Ds q 

The subscript "/i" denotes the two horizontal components of a vector or vectorial operator 
parallel to isobaric surfaces. The quantity ^ = gz is the gravitational potential; k is the 
unit normal to the isobars; / = 2f2 • k is the planetary vorticity, so that fk x vi^ is the 
projection of the Coriolis acceleration onto the isobars; p is the mass density; u = Dp/ Dt is 
the vertical "velocity" in pressure coordinates; s and q are the entropy and radiative heating 
rate per unit mass, respectively; and 



D 

m 



d 



Vh + cj— 

op 



(5) 



is the time derivative following the flow. Equation ([2]) is a statement of hydrostatic equi- 
librium perpendicular to the isobars, so that vertical accelerations have been neglected. 
This is a good approximation on large horizontal scales, since the vertical pressure scale 
height is small compared to the radius of the planet {Hp ~ 10~^R), but it precludes explicit 
representation of small-scale three-dimensional turbulence, which might be an important 
dissipation r nechanism. The continu ity equation (I3l) allows for seepage of material across 
the isobars. IShowman et al.l (l2008bl ) evaluate the heating rate q in the entropy equation 
fror n the vertical divergence q dF/dp of the radiative flux, while IShowman et al.l (l2008al ) 
and iMenou fc Rauscherl (120081 ) use simpler Newtonian schemes in which q oc {T — Teq) / r, 
where Teq is some prescribed reference temperature, which varies with position, and r is a 
prescribed equilibration time. Equation (jl]) is often rewritten in terms of the temperature 
and the specific heat at constant pressure, c^; IShowman et al.l (l2008bl ) rewrite it in terms 
of potential temperature 6 = T{p/poY"'~^'>/"' on the assumption that the adiabatic index 
7 = Cp/cv is constant. What is important here is that q is interpreted to include radiative 
heating only; other sources of entropy such as viscous dissipation have been neglected in the 
simulations that we describe. 



One can derive an expression for the rate of production of mechanical energy based on 



-4- 



the equations above. The first steps are 

^ (1) - ■ ^ - • ^"^ o 

using 1^ . (6) 



L)t p dt 

Recasting uj/p = p^^Dp/Dt as D{pp^^)/ Dt — pD{p^^)/Dt and rearranging terms results in 

Di[r^ + ''^-p)-m=PDi-p- 

The next step is to integrate eq. (JTj) over the mass of the atmosphere down to a depth 
or pressure level sufficiently great that all time derivatives below that level can be neglected, 
but sufficiently shallow so that the vertical gravity g = d^/dz can be treated as spatially 
as well as temporally constant. To do this for the partial time derivative of the potential, 
which is evaluated at constant pressure rather than constant depth, we use eq. ([2]) to write 
the mass element dm = pd^r as dm = g~^dpdA, where dA is an element of area on the 
isobars^ 



\dtjp J \dtjp g dt J g dt J J Dt 

The integral of eq. over the mass of the atmosphere is therefore 

d f f vl p^ 



Using the First Law of Thermodynamics de = Tds —pd{p^^), where e is the internal energy 
per unit mass, we may write eq. fl9ap in the equivalent form 

The integrand on the left side is the Bernoulli constant of the ffow, in which $ does not 
appear because the work done against gravity is balanced by the vertical component of the 
pressure gradient due to the assumption of vertical hydrostatic equilibrium. 



^We have assumed in eq. ([8]) that g'^^dA is constant at a given horizontal position on an isobar as it rises or 
subsides. This would be true in a plane-parallel geometry with strictly constant gravity, since the curvature 
of the isobars would then be second order in the small ratio of the vertical and horizontal lengthscales. But 
in spherical geometry, changes in g and dA can be first order in vertical motions: for example, a spherical 
isobar of varying radius rp{t), has area A (x and gravity gp oc so that g~^dAp (x Vp. We neglect 
this complication since we are mainly concerned with the possibility of steady-state circulation, where the 
isobars would be fixed and {d^/dt)p would vanish a fortiori. 
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3. Prictionless heat engines 

The integrals on the left sides of eqs. iQ are functions of state. Therefore in the absence 
of secular changes, the time averages of the left sides vanish. But the integrands on the 
right sides are not time derivatives of functions of state, and so their averages will not vanish 
automatically. According to equation (jl]), the entropy of each fluid element changes only by 
radiative exchange of heat; generation of heat by viscous dissipation has not been allowed 
for. Insofar as the elements that lose heat are cooler than average (those on the night side), 
and those that gain it are warmer than average (those on the day side), the righthand side 
of equation (l9bl) should be positive. The radiatively driven circulation is a heat engine, 
continuously producing mechanical energy by tapping the flow of heat from the day to the 
night side. 

Real heat engines, whether artificial or natural, have loads: the mechanical energy 
that they produce is ultimately returned to heat through an irreversible process. Engineers 
minimize the dissipation that occurs within the engine itself in order to maximize the power 
transmitted to an external load. The thermodynamic efficiency of the engine is calculated 
as the fraction of the heat absorbed from the hot reservoir that is converted to work rather 
than rejected to the cold reservoir; in this calculation, the heat rejected by the external 
load is not counted. In self-contained natural heat engines, by contrast, dissipation must be 
entirely internal. For example, the kinetic energy of convection is dissipated viscously, even 
at very high Reynolds numbers as in stellar convection zones: turbulent cascades or shocks, 
if convection is transonic, bring kinetic energy down to small scales where viscosity converts 
it to heat. To the extent that a hot Jupiter exchanges energy with its environment purely 
radiatively, so that winds and tides are negligible in the energy balance, and to the extent 
that its energy content is not secularly changing, all of the heat absorbed from the star must 
ultimately be re-radiated. In a sense, therefore, the Carnot efficiency is always zero, but 
only in a sense that would also apply to any heat engine under steady conditions if the heat 
dissipated in the load were counted. 

Equation (jl]) allows mass elements to gain or lose heat by radiation, and as eq. (l9bl) 
shows, these changes can give rise to a net gain in mechanical energy via pdV work. But there 
are no viscous terms in eqs. ([1]) and (j4]) to convert mechanical energy back into heat. There- 
fore, eqs. ([I])-(ll]) describe a frictionless heat engine. The only way that such an engine could 
reach steady state would by adjusting itself to zero Carnot efficiency: that is, by arranging 
for the net pdV work to vanish. An example would be convection so efficient as to erase the 
vertical temperature gradient that drives it. Because the relevant temperature gradient for 
circulation is longitudinal rather than vertical, it is possible that the Carnot efficiency could 
be annulled by a 90° phase shift between the temperature contrast and the radiative heat 
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exchange rather than by complete eradication of longitudinal temperature g radients. Some- 



thing li ke this seems to be going on at the deeper levels of the simulations of IShowman et al. 



( l2008bl ): see their Figure 3. At the shallower levels, however, the temperatures averaged 
over the day and night sides are clearly above and below average, respectively, as might be 
expected from the shorter thermal times there, so a dissipative process is required to explain 
how steady state is reached at these levels. 

Although the dissipation that permits a steady state is not acknowledged in the gov- 
erning equations it may be determined implicitly. If the equations are appropri- 
ately averaged over height /pressure, t hey re d uce t o shallow- water equations [more pre- 
cisely, equivalent-barotropic equations, ISalbyl (Il989l )] with a height-integrated "pressure" 
n oc E*^^'^"^)/'^, where S is the height-integrated density (mass per unit area). It is well 
known that the shallow-water equations permit weak solutions, i.e. nondifferentiable ones, 
that contain dissipative discontinuities ("shocks"); The amount of dissipation in the flow 
as it crosses one of these discontinuities is determined by conservation laws even when the 
dissipative region is unresolved (see, e.g., iLevequd |2002| ). Thus any algorithm that inte- 
grates the equations stably, and that would converge to the correct weak solution (for given 
initial conditions and source terms) in the limit of infinite resolution, will dissipate if the 
weak solution does, eve n if the algorithm doesn't use an explicit dissipative term. Indeed, 
Showman et al.l (l2008bl ) call attention to "hydraulic jumps" in their solutions. 



Unfortunately, the dissipation at hydraulic jumps may be mathematically well-defined 
without being physically correct. That is, the dissipation isn't the same as would be obtained 
from ideal gas dynamics for the same initial conditions and source function, because eqs. ([1])- 
(jl]) are only an approximation to the ideal-gas equations. In particular, the shallow-water 
jump conditions conserve mass and momentum but not energy. The energy per unit mass 
plays a role in the shallow-water approximation that is similar to that pla yed by (—s) in ideal 



gas dynamics: it is c onserved in smooth flow and decreases across shocks (ILandau fc Lifshitz 
19591 : lLevequell2002l ). Inasmuch as eqs. reduce to shallow- water equations when the 

horizontal velocities are constant with height, they probably destroy energy at hydraulic 
jumps. Certainly there is nothing in the equations to convert the dissipated mechanical 
energy to heat. 



4. Other work 



Menou &: Rauscherl (120081 ) adopt equations equivalent to but they add dissipa- 

tive terms. Some of these are diffusive terms for the temperature, vorticity, and velocity- 
divergence fields. The only terms that would appear to dissipate kinetic energy are terms 
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that damp the relative vorticity and velocity divergence; in our notation, these are of the 
forms 

D 

D 

Dt' 

For the purpose of testing their codes against previous calculations made for the Earth's 
atmosphere, they take the friction time Tfric = 27cQ~^, but when they apply the model to 
gaseous exoplanets, they set tmc = oo, perhaps because it is understood that these terms are 
intended to represent turbulent drag against the Earth's solid surface. Ne vertheless, their 



Vh X -t^h 

''"fric 



calcula tions saturate with subsonic wind speeds, whereas the peak speeds in IShowman et al. 



( l2008bl ). and also in most of the previous calculations by Showman and his collaborators, 
are supersonic. 



Showman fc GuillotI (120021 ) integrate the primitive equations in three dimensions. They 
do not include explicit dissipative terms, but they do remark upon a downward energy flux 
at the bottom of their grid, towards the planetary interior. It is not clear whether this is 
enough of an energy sink by itself to explain the quasi-steady circulation in the atmosphere. 



Cho et al.l (120031 ) and lCho et al.l (120081 ) solve two-dimensional shallow- water equations- 



more precisely, equivalent barotropic equations. There is no heat engine in these models: 
thermal "forcing" is represented by a prescribed variation in the effective depth of the cir- 
culating layer, but in such a way that no net work is done on the fluid in the time average. 
For example, the governing equations of I Cho et al.l (120031 ) for the horizontal velocity and 
effective layer thickness, h, are (with some minor changes in notation) 



d 



d 



-gVh{h - ha 
— /C h'Vh • I'h 



+ fkx Vh, 

h-hd 



(CMHS.l) 
(CMHS.2) 



dt I Td 

Here ha and h^ are prescribed functions intended to represent the forcing, and r^, is a constant. 
Integrating the righthand side of eq. (CMHS.l) around a closed streamline yields zero because 
the accelerations parallel to the flow are gradients of scalars. Thus, a steady-solution of these 
equations without dissipative terms could be compatible with energy conservation|§ In fact. 



^This argument docs not apply to ([T]) because streamlines may cross isobaric surfaces. Thus if ABB' A' is 
a hypothetical "rectangular" streamline such that the legs AB and B'A' lie entirely on neighboring isobars 
p and p', then the line integral of the three-dimensional gradient V<i> does vanish, but the contributions of 
the vertical legs BB' and A' A are offset by the pressure acceleration due to hydrostatic equilibrium [eq. 
and these contributions do not cancel unless pbb' = PA' A- 
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a steady solution with = and h = ha clearly exists if hd = ha- On the other hand, the 
equations do not conserve the mass in the layer unless —>■ oo, in which case (CMHS.2) can 
be recast as dTi/dt + Vh • (Svh) = for areal mass density S oc h^^'^. 



Langton fc LaughlinI (120071 ) integrate the shallow-water equations using a spectral method. 



They do not display any explicit dissipative term in th e morn entum equations, but they men- 
tion the use of a hyperviscosity. iLangton fc LaughlinI (120081 ) adopt somewhat different, but 
also two-dimensional, equations with a radiative forcing that appears to yield a nonconser- 
vative force in the momentum equation [the term i?T Vh In p in their eq. (9)], so they may 
have a heat engine. Their only explicit mechanical dissipation is again a hyperviscosity, 
—BV^V\i, but the mechanical energy lost via this term is apparently not added as heat to 
the evolutionary equation for the temperature. 

Not all of the work in this subject ha s been carried out wit h th e primitive, shallow-water, 
or equivalent-barotropic approximations. iBurkert et al.l (120051 ) and lDobbs-Dixon fc LinI (120081 ) 
use the full ideal-gas equations, including radiative forcing terms but without the approx- 
imation of vertical hydrostatic equilibrium. The former work is limited to two dimensions 
(radius and longitude), while the latter is fully three dimensional. It is not clear from their 
abbreviated descriptions whether their schemes conserve total energy, but if they do, then 
these schemes should be correctly converting the mechanical energy of the circulation into 
heat. Even so, by n eglect ing explicit dissipation, they may overestimate the flow speed. 
Dobbs-Dixon &: LinI (120081 ) remark upon flow speeds up to Mach 2.7 and dissipation in 
shocks. 



5. Discussion 

Analogies between atmospheric circulation and heat engines have been made before. 
In the atmospheres of planets such as the Earth whose rotational periods are short com- 
pared both to their orbital periods and to the thermal times in their atmosphere^, the 
primary temperature gradient is latitudinal, and the atmospheric heat engine operates be- 
tw een the equator and the poles, via the Hadley, mid-latitude, and polar cells. As discussed 



by iLorenz et al.l (120011 . and references therein), the latitudinal circulation on Earth, and 



perhaps also on Mars and Titan, adjusts itself approximately so as to maximize the rate of 
entropy production, defined as 5 = (5adv(^poio ~ ^equator)) where Qadv is the rate of poleward 
heat transport by advection. Here 5* does not represent a secular increase in the entropy of 



•^The temperature of the terrestrial atmosphere varies rather httle over the diurnal cycle except at low 
altitude where the air easily exchanges heat with the ground. 



- 9 - 



the atmosphere itself, since the heat and entropy brought to the poles is ultimately radiated 
to space, but rather a (small) increase in the entropy of the universe. Maximizing S is not 
equivalent to maximizing Carnot efficiency (?7c); but neither is it compatible with rjc = 0. 
The entropy production defined this way would vanish for Qadv = 0, which would lead to 
the largest difference between equatorial and polar temperatures and hence to the largest 
possible rjc] but it would also vanish if the heat transport were so large that Tpoie = Tequator, 
when ?7c 0. In the actual state of approximately maximal S, the mechanical energy of the 
circulation is continuously produced by the heat engine but dissipated by turbulent friction 
with the solid surface. Since the transported heat is ultimately radiated, a quantity related 
to S that could be measured in those simulations that include radiative transport might be 

photosph. 

where -Frad is the difference between the outward planetary flux and the inward stellar 
flux, and Tph is the gas temperature at the photosphere. Since the optical and infrared 
photospheres do not coincide, however, a better definition might be 



'LL — 

atmos. 




- rad 



dV, (10) 



where the subscript indicates that this is entropy production in the sense of iLorenz et al. 



(I2OOII ). Normalized perhaps by the total insolation and by a mean photospheric temperature, 
the quantity (fTOl) might be used to assess the degree to which the redistribution of heat is 
maximal. It is not obvious a priori why S'll should be maximal, however. 

There are at least two distinct though related ways in which neglect of explicit dissipation 
could lead to erroneous results: 



1. If mechanical energy is dissipated by numerical effects that do not convert it to heat, 
then energy is lost. The seriousness of this error for predictions of infrared light curves, 
etc., depends upon the fraction of the advected energy that is carried in kinetic form 
rather than as enthalpy. The latter could be measured by integrating the two contribu- 
tions to the energy fiux over the meridional plane at the dawn and dusk terminators, 
for example. 

2. Even when total energy is conserved or its loss is unimportant, wind speeds and rates 
of redistribution of heat may be overestimated when physical mechanisms of drag on 
the fiow are left out. 
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Terrestrial circulation is highly subsonic, so modelers presume that errors of the first type 
are not serious. This is not so obvious for the models of hot Jupiters. Since the simulators 
do not report the balance between the kinetic and enthalpic energy fluxes, consider the 
hypothesis that the wind speed U at each (sufficiently deep) pressure level adjusts itself so 
that the advection time Tadv = R/U is proportional to the thermal time at that depth, 

_ p npHp KP^Hp 
Tth = c " gaT^ ' 

where Hp ^ k-oT /2mpg is the pressure scale height. In an isothermal atmosphere, this 
hypothesis predicts U oc k~ V^^. The horizontal enthalpy flux would then scale with depth as 
and the kinetic flux as i^~^p~^ . Since the opacity generally increases with pressure, the 
depth-integrated heat transport would be dominated by low pressures near the photosphere, 
where the kinetic component is important inasmuch as the speeds are transonic according 
to the simulations. 

On the other hand, if the simulations have overestimated wind speeds, then the kinetic 
contribution to the energy flux may turn out to be minor after all. This brings us to 
point 2 above: What physical mechanisms of dissipation and drag may operate in real 
exoplanetary atmospheres that have not been included in the simulations? Unlike terrestrial 
planets, hot Jupiters presumably have no solid surface to provide boundary-layer drag; a 
dissipative mechanism that is distributed through the atmosphere is required. Two important 
possibilities are shocks and shear-driven turbulence (perhaps also waves). Shocks would 
occur in steady flow only if it is supersonic. They might do little to mix the atmosphere 
and to resist chemical fractionation. If shocks dominate, this woul d be good for s i mulat ions, 



which ought to be able to capture them, as in the simulations of iBurkert et al.l (120051 ) and 



Dobbs-Dixon fc Lid (120081 ). If turbulence dominates, however, then the simulations may 
require a subgrid model, because a necessary characteristic of the turbulence is that it extend 
to unresolvably small scales where viscosity is effective. The subgrid model should be based 
on an understanding of the instabilities that promote the turbulence and how they saturate, 
unless an empirical model can be calibrated against the well-resolved planetary atmospheres 
of the solar system. 

The principal obstacle to the development of three-dimensional turbulence and a cascade 
to small scales is the stable stratification. The kinetic energy available in the vertical wind 
shear must be greater than the potential energy required to mix the atmosphere. The ratio 
of these energies is quantified by the Richardson number Ri = N'^ / {dU / dr)'^ , where is the 
Brunt-Vaisala frequency and U is the horizontal velocity. A requirement for KH instability 
in vertical planes i s that Ri < 1/4 som ewhere in the flow, at least under adiabatic and 



inviscid conditions (jPrazin fc Reidlll98ll ). It would be helpful if numerical researchers would 
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report the values of Ri achieved in their simulations. For an isothermal atmosphere with 
density and pressure scale height Hp = c'^/'yg, the condition Ri < ^ is equivalent to 



dU 



dr 



> 



1 



so with 7 ^ 7/5 as for nondegenerate molecular hydrogen, the flow must be transonic or else 
change on a scale smaller than the scale height in order that Ri < 1/4. But it is possible 
that instability may occur at larger Richardson number in the presence of radiative dif f usion. 



Othe r instabilities that may be relevant are the b aroclinic instabilitv (ll^edlosky 



20061 ) ■ the Goldreich-Schubert-Fricke instability (IGoldreich fc SchubertI 



and perhaps even the magnetorotational instability (iBalbus &: Hawley 



967 



1987 



Fricke 



Vallis 



19681 



19911 ) if the shear 



extends to such a depth that the atmosphere becomes substantially conducting, and if the 
shear has the right sign. It is beyond the scope of this note to assess the importance of these 
instabilities for controlling the speed of the circulation. Nevertheless, we will discuss the 
GSF instability briefly because it is illustrative and must surely occur at some level. 

The GSF instability is enabled by thermal diffusion in baroclinic atmospheres of neg- 
ligible viscosity that do not rotate on cylinders, d Q/dz ?^ 0, ev en if the atmosphere is 
dynamically stable according to H0lland's Criterion (lTassoullll978l ). i.e. baroclinically sta- 
ble. The instability is axisymmetric, and the maximum growth rate is expressed in terms of 
the cylindrical radius cc and angular momentum per unit mass j = w'^Q by 



2^3 



dw 



1/2 



but this is approached only at wavelengths A < 2'k{x/ ^Y^"^ ■, where x = 16(tT^/3kp^c„ ~ 
{c/ Kp){pT:ad/p) is the thermal diffusivity, so that radiative diffusion undercuts the restoring 
force of buoyancy. Thus A < 1 km (p/bar)~^(K/cm^ g^^)^^/^(T/10^K)^. Such a wavelength 
would not be resolved by global simulations. One might expect the GSF instabilities to 
saturate nonlinearly at displacements ^ ~ A due to nonaxisymmetric Kelvin-Helmholtz in- 
stabilities acting on the rising and fallin g "fingers," which carry oppo sitely signed eulerian 
angular-velocity perturbations. However, IGoldreich &: SchubertI (119671 ) themselves took the 
point of view that the saturation occurs with displacements comparable to the pressure 
scale height, and therefore provides a turbulent viscosity ~ QHp/27c. This is something 
best studied by local rather t han g lobal simulations, and then represented in the latter by a 
subgrid model. iKorycanskyI (1199 ll ) simulated the nonlinear outcome of the GSF instability 



in the special case where dj /dz = and A^^ > —zu ^dj'^/dw > 0. His results support 



Goldreich fc SchubertI (Il967l )'s view of the saturation, but his two-dimensional simulations 
could not have represented the nonaxisymmetric Kelvin-Helmholtz instabilities most likely 
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to limit GSF modes. lArlt &: UrpinI (120041 ) simulated the case dy^ jdw > and dj'^/dz < in 
three dimensions, using ZEUS3D, and concluded that mixing was efficient, but they adopted 
a bar otropic equat i on of state so that their unperturbed state had to be out of equilib- 



Menou et al.l (120041 ) added magnetic effects to the GSF analysis, but only in the linear 



rmm. 
regime. 

While the growth rate of instabilities that rely on thermal diffusion should decrease 
rapidly with increasing pressure, the rate at which radiative transfer tends to restore the 
stratification also decreases, so that the outcome for the proffies of entropy and angular ve- 
locity is unclear. If any of these instabilities is effective at redistributing angular momentum, 
then the radiatively driven circulation may go deeper than present simulations suggest, and 
the time required for the rotation proffie to reach steady state may be very long. 



We thank Adam Burrows, Kristen Menou, Jonathan Mitchell, Geoffrey Vallis, and the 
Peyton-Hall astro-ph coffee klatsch for helpful discussions. This work was supported in part 
by the National Science foundation under grant AST-0707373. 
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